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Abstract. We study a class of optical circuits with vacuum input states 
consisting of Gaussian sources without coherent displacements such as down- 
converters and squeezers, together with detectors and passive interferometry 
(beam-splitters, polarisation rotations, phase-shifters etc.). We show that the 
outgoing state leaving the optical circuit can be expressed in terms of so-called 
multi-dimensional Hermite polynomials and give their recursion and orthogonality 
relations. We show how quantum teleportation of photon polarisation can be 
modelled using this description. 

Suppose we have an optical circuit, that is, a coUection of various optical 
components. It is usually important to know what the outgoing state of this circuit is. 
In this paper, we give a description of the outgoing state for a special class of optical 
circuits with a special class of input states. 

First, in section |l], we define this class of optical circuits and show that the; 
can be described by so-called multi-dimensional Hermite polynomials. In section ^ 
we give an example of this description. Section ^ discusses the Hermite polynomials 
and finally, in section ^, we briefly consider the effect of imperfect detectors on the 
outgoing state. 

1. The Optical Circuit 

What do we mean by an optical circuit? We can think of a black box with incoming 
and outgoing modes of the electro-magnetic field. The black box transforms a state 
of the incoming modes into a (different) state of the outgoing modes. The black box 
is what we call an optical circuit. We can now take a more detailed look inside the 
black box. We will consider three types of components. 

First, the modes might be mixed by beam-splitters, or they may pick up a 
relative phase shift or polarisation rotation. These operations all belong to a class of 
optical components which preserve the photon number. We call them passive optical 
components. 

Secondly, we may find optical components such as lasers, down-converters or 
(optical) parametric amplifiers in the black box. These components can be viewed as 
photon sources, since they do not leave the photon number invariant. We will call 
these components active optical components. 

And finally, the box will generally include measurement devices, the outcomes 
of which may modify optical components on the remaining modes depending on 
the detection outcomes. This is called feed-forward detection. We can immediately 
simplify optical circuits using feed-forward detection, by considering the family of 
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fixed circuits corresponding to the set of measurement outcomes (see also Ref. Q). 
In addition, we can postpone the measurement to the end, where all the optical 
components have 'acted' on the modes. 

These three component types have their own characteristic mathematical 
description. A passive component yields a unitary evolution Ui , which can be written 
as 

= cxp -IK ^ Cjkajal - H.c. , (1) 

V 3k J 

where H.c. denotes the Hermitian conjugate. This unitary evolution commutes with 
the total number operator n = o.'jO.j- 

Active components also correspond to unitary transformations, which can be 
written as exp{—itHj-'^). Here i/j-'' is the interaction Hamiltonian associated with the 
j^^ active component in a sequence. This Hamiltonian does not necessarily commute 
with the total number operator. To make a typographical distinction between passive 
and active components, we denote the i^^ passive component by Ui, and the j'^ active 
component by its evolution in terms of the interaction Hamiltonian. 

The mathematical description of the (ideal) measurement will correspond to 
taking the inner product of the outgoing state prior to the measurement with the 
eigenstate corresponding to the measurement. 

Now that we have the components of an optical circuit of N modes, we have to 
combine them into an actual circuit. Mathematically, this corresponds to applying 
the unitary evolutions of the successive components to the input state. Let IV'in) be 
the input state and IV'prior) the output state prior to the measurement. We then have 
(with K > some integer) 



l^prior) = f/i^e-'*^^ ...C/ie-'*«^ f/olV'in), (2) 

where it should be noted that Ui might be the identity operator 1 or a product of 
unitary transformations corresponding to passive components: 

U^^ll U,,k . (3) 
k 

When the (multi-mode) eigenstate corresponding to the measurement outcome for a 
limited set of modes labelled 1, . . . , M with M < TV is given by I7) = |ni, 712, • ■ • , "a/) 
with M the number of detected modes out of a total of N modes, and rii the number 
of photons found in mode i, the state leaving the optical circuit in the undetected 
modes is given by 

\^out)M+l,...,N =1,...,M (7lV'prior)l,...,Ar . (4) 

In this paper, we study the outgoing states IV'out) for a special class of optical 
circuits. First, we assume that the input state is the vacuum on all modes. Thus, we 
effectively study optical circuits as state preparation devices. Secondly, our class of 
optical circuits include all possible passive components, but only active components 
with quadratic interaction Hamiltonians: 

H^/'=Y^alR^^,^aj+Y^a,R^td,, (5) 

kl kl 

where i?*-^^ is some complex symmetric matrix. This matrix determines the behaviour 
of the j^^ active component, which can be any combination of down-converters 
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and squeezers. Finally, we consider ideal photo-detection, where the eigenstate 
corresponding to the measurement outcome can be written as I7) = |ni, . . . jUm)- 

The class of optical circuits we consider here is not the most general class, but 
it still includes important experiments like quantum teleportation |Q, entanglement 
swapping and the demonstration of GHZ correlations Q . In section ^ we show how 
teleportation can be modelled using the methods presented here. 

The state jV') prior to the photo-detection can be written in terms of the 
components of the optical circuit as 

1^) = Ukb-'^^'''''' . . . C/ie"**^^" |0) . (6) 

The creation and annihilation operators and for mode i satisfy the standard 
canonical commutation relations 

[flj, a]] = 6ij and [ai,dj] = [a]', a]] = , (7) 

with j = 1 . . . N. 



For any unitary evolution U , we have the relation 



E 



II 



E 

(=0 



{uRwy 
II 



„URW 



(8) 



Furthermore, if U is due to a collection of only passive components, such an evolution 
leaves the vacuum invariant: [/|0) = |0). Using these two properties it can be shown 
that Eq. (0) can be written as 



I'/') = exp 



1 



N 



2 ^ 



T 



1 ^ 

- ^ OiAijCij 



|0) 



(9) 



where A is some complex symmetric matrix. We will now simplify this expression by 
normal-ordering this evolution. 

Define {a,Aa) = diAijdj. As shown by Braunstein [^, we can rewrite Eq. 
using two passive unitary transformations U and V as: 

where A is a diagonal matrix with real non-negative eigenvalues A;. This means that, 
starting from vacuum, the class of optical circuits we consider here is equivalent to 
a set of single-mode squeezers, followed by a passive unitary transformation U and 
photo-detection. Since A is diagonal, we can write Eq. (10) as 

/ N 



|0> 



(11) 



We can now determine the normal ordering of every factor exp[— ^(a|^)^ 4- ^df] 
separately. Note that the operators af , (a|^)^ and 2a|ai-t-l generate an su(l,l) algebra. 
According to Refs. P, Q, H], this may be normal-ordered as 



g — Ai tanh I -J- 1 (dj) g — 2 ln(cosh | -^^ | )d] d^ ^A^ tanh | -^^ |dj (12) 



where = Ai/|Ai|. In general, when L± and Lq are generators of an su(l, 1) algebra 
(i.e., when A is unitary) we find S 



i(TL+r-L_) 



r tanh \t\Lj^ ^ — 2 In (cosh \ t\)Lq tanh |t|L_ 



(13) 
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with r a complex couphng constant and f its orientation in the complex plane. When 
we now apply this operator to the vacuum, the annihilation operators will vanish, 
leaving only the exponential function of the creation operators. We thus have 

^ C/e-^('^*'^*^*V^|0) = e-5(snB'J+)|o) , (14) 

with B = UA*U'' , again by virtue of the invariance property of the vacuum. This is 
the state of the interferometer prior to photo-detection. It corresponds to multi-mode 
squeezed vacuum. 

The photo-detection itself can be modelled by successive application of 
annihilation operators. Every annihilation operator removes a photon in mode 
i from the state Suppose the optical circuit employs N distinct modes. We 

will now detect M modes, finding ni + . . . + um = ^tot photons (with M < N). 
These modes can be relabelled 1 to Af. The vector n denotes the particular detector 
'signature': n = (rii, . . . ,nA/) means that ni photons are detected in mode 1, n2 in 
mode 2, and so on. The freely propagating outgoing state \ipfi) can then be described 
as 

IV'ri) = i..M(ni, • • ■ , 7^A^|■0)l..JV = c,-j(0|a"' • • • a^f 



(15) 



Here, cn = {nil ■ ■ ■ um^) ^ . 

At this point we find it convenient to introduce the A'^-mode Bargmann 
representation The creation and annihilation operators obey the commutation 
relations given in Eq. (|^). We can replace these operators with c- numbers and their 
derivatives according to 

' (16) 



ai 



dai 



and Oi — > di 

The commutation relations then read 

[di, aj] = Sij and [di, dj] = [a^, a-j] = . (17) 

Note that the actual values of ai are irrelevant (the creation and annihilation operators 
do not have numerical values either) ; what matters here is the functional relationship 
between ai and 9q . . 

The state created by the optical circuit in this representation (prior to the 
detections, analogous to Eq. (14)) in the Bargmann representation is 



V'(cf) = exp 



-[a,Ba) 



exp 



1 ^ 

-^a,B,jaj 

ij 



(18) 



Returning to Eq. ( [T5| ) , we can write the freely propagating state after detection of the 
auxiliary modes in the Bargmann representation as 



Md) ^ an d^^ ■ ■ ■ dlr 



g-i(a,Ba) 



(19) 



up to some normalisation factor, where d' = {ai, . . . ,aM)- By setting o?' = we 
ensure that no more that Ui photons are present in mode i. It plays the role of the 
vacuum bra in Eq. (p^). 

Now that we have an expression for the freely propagating state emerging from 
our optical setup after detection, we seek to simplify it. We can multiply ipfiid) by 
the identity operator I, written as 



/=(-l) 



2A't. 



■ exp 



" 1 / - ^; 




'1/^ " 


--{a,Ba) 


exp 


-[a,Ba} 



(20) 
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where TVtot is the total number of detected photons. We then find the foUowing 
expression for the unnormahsed freely propagating state created by our optical circuit: 

M^) « Cn{-lf'°'H?{a) e-^("^^'S) . (21) 

q'— 

Now we introduce the so-called multi- dimensional Heirmte polynomial, or MDHP 
for short: 

The use of multi-dimensional Hermite polynomials and Hermite polynomials of two 
variables have previously been used to describe A^-dimensional first-order systems 
To[ |ll] and photon statistics |12, ^ Here, we have shown that the lowest order 



of the outgoing state of optical circuits with quadratic components (as described by 
Eq. (|^)) and conditional photo-detection can be expressed directly in terms of an 
MDHP. 

In physical systems, the coupling constants (the A^'s) are usually very small 
(i.e., Ai ^ 1 or possibly Ai < 1). This means that for all practical purposes 
only the first order term in Eq. ( pT] ) is important (i.e., for small A^'s we can 
approximate the exponential by 1). Consequently, studying the multi-dimensional 
Hermite polynomials yield knowledge about the typical states we can produce using 
Gaussian sources without coherent displacements. In section ^ we take a closer look 
at these polynomials, but first we consider the description of quantum teleportation 
in this representation. 



2. Example: Quantum Teleportation 

As an example of how to determine the outgoing state of an optical circuit, consider the 
teleportation experiment by Bouwmeester et al. |^ . The optical circuit corresponding 
to this experiment consists of eight incoming modes, all in the vacuum state. 
Physically, there are four spatial modes a, b, c and d, all with two polarisation 
components x and y. Two down-converters create entangled polarisation states; they 
belong to the class of active Gaussian components without coherent displacements. 
Mode a undergoes a polarisation rotation over an angle 9 and modes b and c are 
mixed in a 50:50 beam-splitter. Finally, modes b and c emerging from the beam- 
splitter are detected with polarisation insensitive detectors and mode a is detected 
using a polarisation sensitive detector. The state which is to be teleported is therefore 
given by 

1^-) = cos 6l|a;) -sin 6*12/) . (23) 

The state prior to the detection and normal ordering (corresponding to Eq. (||)) 
is given by (r is a coupling constant) 

with 



V2 



/ 1 \ 

10 

10 

\ 1 / 



(25) 
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and u'f = (al.,a^,bl.,b^y), = (c|., , rft ^ rft ). This can be written as 



^{d\Ad^) + ^{d, Ad) 



with d = {dx , ■ 



l^prior) = exp 
■ ,dy) and A the (symmetric) matrix 



|0) 



(26) 



A = 



V2 



/ 





— sin cos 9 


— sin 9 


cos 9 





\ 







cos 6 sin 9 


cos 9 


sin 9 
























-1 















-1 




















1 













1 





















V 












/ 



(27) 



We now have to find the normal ordering of Eq. (|2q). Since A is unitary, the 
polynomial (a^, Ad^) is a generator of an su(l, 1) algebra. According to Truax the 
normal ordering of the exponential thus yields a state 



l^pri. 



exp 



-(a^^at) 



|0) 



(28) 



with ^ — (r tanh |t|)/|t|. The lowest order contribution after three detected photons 
is due to the term ^^(a^, yl.a^^)^/8. However, first we write Eq. ( p8| ) in the Bargmann 
representation: 



V'prior(a) 



exp 



•-{a, Ad) 



(29) 



where cJ = (aa^ , . . . , ) and cJ' = (aa^ , . . . , Ofc ) . The polarisation independent 
photo-detection is then modelled by the differentiation {db^dcy — db^dcj)- Given a 
detector hit in mode a^, the polarisation sensitive detection of mode a is modelled by 



V'out((S) = da^ {db^dcy - db^dc^) exp 



-(d, Acx) 



S'=0 

The outgoing state in the Bargmann representation is thus given by 



■0out (a) 



I COS f aA 



sm U ad^ 



(30) 
(31) 



which is the state teleported from mode a to mode d in the Bargmann representation. 
This procedure essentially amounts to evaluating the multi-dimensional Hermite 
polynomial H^{d). Note that the polarisation independent detection of modes b 
and c yield a superposition of the MDHP's. 



3. The Hermite Polynomials 

The one-dimensional Hermite polynomials are of course well known from the 
description of the linear harmonic oscillator in quantum mechanics. These polynomials 
may be obtained from a generating function G. Furthermore, there exist two recursion 
relations and an orthogonality relation between them. The theory of multi-dimensional 
Hermite polynomials with real variables has been developed by Appell and Kempe de 
Feriet [ p^ and in the Bateman project |T^. Mizrahi derived an expression for real 
MDHP's from an n-dimensional generalisation of the Rodriguez formula We will 
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now give the generating function for the complex MDHP's given by Eq. ( p^ ) and 
consecutively derive the recursion relations and the orthogonality relation (see also 
Ref. 

Define the generating function GB{a,f3) to be 

Gsid, P) = e("^^«-^(^~.^« = y ^ . . . ^Hi{d) . (32) 

n 

Gsia, /?) gives rise to the MDHP in Eq. (|2^), which determines this particular choice. 
Note that the inner product (a, BP) does not involve any complex conjugation. If 
complex conjugation was involved, we would have obtained different polynomials 
(which we could also have called multi-dimensional Hermite polynomials, but they 
would not bear the same relationship to optical circuits). 

In the rest of the paper we use the following notation: by n — Cj we mean that the 
j**^ entry of the vector n = (ni, . . . , um) is lowered by one, thus becoming rij — 1. By 
differentiation of both sides of the generating function in Eq. ( p2| ) we can thus show 
that the first recursion relation becomes 

—Hi id) = J2 B^,n,Hi_,^ (a) . (33) 
The second recursion relation is given by 

M M 

H,i+e. («) - E + E B^,n,H?_^^ (a) = , (34) 

which can be proved by mathematical induction using 

M M 

Y,B,knkH?_,^^,^{a) - BuH?(d) = Y,B^kmkH?+,^{d) . (35) 
fc=i fe=i 

Here, we have chosen m — n — Ck- 

The orthogonality relation is somewhat more involved. Ultimately, we want to 
use this relation to determine the normalisation constant of the states given by Eq. 
(^T|). To find this normalisation we have to evaluate the integral 

ddil;^{d)tp^{d) . (36) 

The state V'n includes |<;'=o, which translates into a delta- function S{d') in the 
integrand. The relevant integral thus becomes 

/ da e-^^("'^") [H?id)] * H?id) <5(d') . (37) 

From the orthonormality of different quantum states we know that this integral must 
be proportional to (Jfj.^. 

Since in the Bargmann representation we are only concerned with the functional 
relationship between and da^ and not the actual values, we can choose to be 
real. To stress this, we write ai Xi. The orthogonality relation is thus derived from 

dxr^mraix)^ [ dxe~('^'''<^^'^HS'{x)Hg{x)5{x') . (38) 
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where 5{x') is the real version of S{a'). FoUowing Klauderer [|ll| we find that 



(39) 



where is the differential operator dJ^^ ■ ■ ■ d"^l acting solely on the exponential 
function. We now integrate the right-hand side by parts, yielding 



/+00 

_(_l)JVt„t /'dfe-3(^'^^)a|-"'e-^(^'^*^)a,.i75(f) , (40) 



Xi — — 00 



with d'x — dxi ■ ■ ■ dxi^idxi^i ■ ■ ■ dxpf. The left-hand term is equal to zero when 
Re(i?) is positive definite, i.e., when (x,Ke{B)x) > for all non-zero x. Repeating 
this procedure rii times yields 



e 2 



iS,B'x)Qn.^jjB^^^ . (41) 



When there is at least one > mi , differentiating the MDHP rii times to Xi will yield 
zero. Thus we have 

J dxe-'-^'^'^^^^^'^Hf {x)Hg{x) = for n ^ m (42) 

when Re(-B) is positive definite and rii ^ nii for any i. The case where n equals m is 
given by 

j e-5(-~ R<=(^)-)iJ?* {x)Hg{x) = SnrnAf , (43) 

where Sfi^ denotes the product of Smrm with 1 < * < Here, Af is equal to 

AA= 2^-'B;'ii---B]^^ni!---nAr! |7r-iBp^ . (44) 
For the proof of this identity we refer to Ref. ||ll| . 



4. Imperfect Detectors 

So far, we only considered the use of ideal photo-detection. That is, we assumed that 
the detectors tell us exactly and with unit efficiency how many photons were present 
in the detected mode. However, in reality such detectors do not exist. In particular 
we have to incorporate losses (non-perfect efficiency) and dark counts. Furthermore, 
we have to take into account the fact that most detectors do not have a single-photon 
resolution (i.e., they cannot distinguish a single photon from two photons) [Q. 

This model is not suitable when we want to include dark counts. These unwanted 
light sources provide thermal light, which is not of the form of Eq. (^). In single-shot 
experiments, however, dark counts can be neglected when the detectors operate only 
within a narrow time interval. 
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We can model the efRciency of a detector by placing a beam splitter with 
transmission amplitude 77 in front of a perfect detector |Q . The part of the signal 
which is reflected by the beam-splitter (and which will therefore never reach the 
detector) is the loss due to the imperfect detector. Since beam splitters are part 
of the set of optical devices we allow, we can make this generalisation without any 
problem. We now trace out all the reflected modes (they are truly 'lost'), and end up 
with a mixture in the remaining undetected modes. 

Next, we can model the lack of single-photon resolution by using the relative 
probabilities p{n\k) and p{m\k) of the actual number n or to of detected photons 
conditioned on the indication of k photons in the detector (as described in Ref. ) . 
We can determine the pure states according to n and m detected photons, and add 
them with relative weights p{n\k) and p{m\k). This method is trivially generalised for 
more than two possible detected photon numbers. 

Finally, we should note that our description of this class of optical circuits (in 
terms of multi-dimensional Hermite polynomials) is essentially a one-way function. 
Given a certain setup, it is relatively straightforward to determine the outgoing state 
of the circuit. The other way around, however, is very difficult. As exemplified by 
our efforts in Ref. [Q , it is almost impossible to obtain the matrix B associated with 
an optical circuit which produces a particular predetermined state from a Gaussian 
source. 



5. Conclusions 



In this paper, we have derived the general form of squeezed multi-mode vacuum states 
conditioned on photo-detection of some of the modes. To lowest order, the outgoing 
states in the Bargmann representation are proportional to multi-dimensional Hermite 
polynomials. As an example, we showed how teleportation can be described this way. 
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